A&A manuscript no. 

(will be inserted by hand later) 



Your thesaurus codes are: 

08.01.1, 08.01.3, 08.08.2, 08.09.2, 08.15.1, 08.18.1 



ASTRONOMY 
AND 
ASTROPHYSICS 
1.2.2008 



OS 
OS 

os 



Spectral analysis of the multi mode pulsating subluminous 
B star PG 1605+072 

U. Heber 1 , I.N. Reid 2 , K. Werner 3 

1 Dr. Remeis-Sternwarte, Astronomisches Institut der Universitat Erlangen-Niirnberg, Sternwartstr. 7, D 96049 Bamberg, Ger- 
many (E-mail heber@sternwarte.uni-erlangen.de) 

2 Palomar Observatory 105-24, California Institute of Technology, Pasadena, CA 91125, USA 

3 Institut fur Astronomie und Astrophysik, Universitat Tubingen, D 72076 Tubingen 

Received ; accepted 



o 

(N 



> 

<n : 

SO . 
(N ■ 

■ 

o : 

Os . 
Os ■ 

^ ' 

Or 

6 ■ 
a: 



Abstract. PG 1605+072 has unique pulsational proper- 
ties amongst the members of the new class of pulsat- 
ing sdB (EC 14026) stars. It has the longest periods and 
the richest, most puzzling frequency spectrum (55 peri- 
ods) . We present a quantitative analysis of a Keck HIRES 
spectrum using NLTE and LTE model atmospheres. At- 
mospheric parameters ( T c g, log g, log^s.), metal abun- 
dances, and the rotational velocity are determined. He, C, 
N and O are subsolar as well as the intermediate mass 
elements Mg and Si, whereas Ne and Fe are solar. This 
abundance pattern is caused by diffusion. PG 1605+072 
displays considerable line broadening probably caused by 
rotation (vsini = 39km/s, P<8.7h) and, therefore, is pre- 
dicted to evolve into an unusually fast rotating white 
dwarf. Unequal rotational splitting may explain its puz- 
zling pulsation pattern. The solar Fe abundance and the 
rapid rotation nicely confirms recent predictions of diffu- 
sion/pulsation models. 

Key words: stars: abundances — stars: atmospheres 
— stars: oscillations — stars: horizontal branch — stars: 
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1. Introduction 

Hot subluminous B stars (sdB) form a homogeneous group 
dominating the population of faint blue stars (B<16 m ). 
Following ideas outlined by Heber (1986) the sdB stars 
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can be identified with models for extreme HB (EHB) stars, 
which differ markedly from those for normal HB stars. 
An EHB star bears great resemblance to a helium main- 
sequence star of half a solar mass and its further evolution 
should proceed similarly (i.e. directly to the white dwarf 
graveyard, Dorman et al. 1993). 

Recently, several sdB stars have been found to be pul- 
sating (termed EC14026 stars after the prototype, see 
O'Donoghue et al. 1999 for a review), defining a new in- 
stability strip in the HR-diagram. The study of these pul- 
sators offers the possibility of applying the tools of as- 
teroseismology to investigate the structure of sdB stars. 
The existence of pulsating sdB stars was predicted by 
Charpinet et al. (1996), who uncovered an efficient driving 
mechanism due to an opacity bump associated with iron 
ionization in EHB models. However, in order to drive the 
pulsations, iron needed to be enhanced in the appropriate 
subphotospheric layers, possibly due to diffusion. Subse- 
quently, Charpinet et al. (1997) confirmed this assumption 
by detailed diffusion calculations. Even more encouraging 
was the agreement of the observed and predicted instabil- 
ity strip. 

Thirteen pulsating sdB stars are well-studied photo- 
metrically (O'Donoghue et al. 1999). A precise knowledge 
of effective temperature, gravity, element abundances and 
rotation is a prerequisite for the asteroseismological inves- 
tigation. 

PG 1605+072 was selected for a detailed quantitative 
spectral analysis because it displays the richest frequency 
spectrum amongst the EC 14026 stars (>50 periods have 
been identified, Kilkenny al. 1999). Recently, Kawaler 
(1999) predicted from his modelling of the pulsations that 
PG 1605+072 should be rotating. 

2. Observation and data reduction 

A high resolution optical spectrum of PG 1605+072 was 
obtained with the HIRES echelle spectrograph (Vogt et 
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al. 1994) on the Keck I telescope on July 20, 1998 using 
the blue cross disperser to cover the full wavelength region 
between 3700A and 5200A at a resolution 0.09A. 

The exposure time (600s) is longer than the pulsational 
periods (206-573s). The standard data reduction as de- 
scribed by Zuckerman & Reid (1998) resulted in spectral 
orders that have a somewhat wavy continuum. In order to 
remove the waviness we used the spectrum of H1504+65 (a 
very hot pre-white dwarf devoid of hydrogen and helium, 
Werner 1991) which was observed in the same night. Its 
spectrum has only few weak lines of highly ionized met- 
als in the blue (3600-4480A) where the strong Balmcr 
lines are found in the sdB stars. Therefore we normal- 
ized individual spectral orders 1 to 20 (3600-4480A) of 
the sdB stars by dividing through the smoothed spectrum 
of H1504+65. The remaining orders were normalized by 
fitting the continuum with spline functions (interpolated 
for orders 26 and 27 which contain H/3). Judged from the 
match of line profiles in the overlapping parts of neigh- 
boring orders this procedure worked extremely well. At- 
mospheric parameters determined from individual Balmer 
lines are found to be consistent with each other except for 
H/3. Therefore, we excluded H/3 from the fit procedure. 
Moreover, the resulting T c g and log g are also in excel- 
lent agreement with those from the fit of a low resolu- 
tion spectrum obtained at the ESO NTT (provided by S. 
Moehler). Weak lines of C, N, O, Ne, Mg, Si, and Fe can 
also be identified. 

3. Atmospheric parameters 

The simultaneous fitting of Balmcr and He line profiles 
by a grid of synthetic spectra (see Saffer et al. 1994) has 
become the standard technique to determine the atmo- 
spheric parameters of sdB stars. The Balmer lines (H7 to 
H12), He 1 (4471 A, 4026A, 4922A, 4713A, 5016A, 5048A) 
and He 11 4686A lines are fitted to derive all three param- 
eters simultaneously. 

The analysis is based on grids of metal line blanketed 
LTE model atmospheres for solar metalicity and Kurucz' 
ATLAS6 Opacity Distribution Functions (see Heber et al. 
1999). Synthetic spectra are calculated with Lemke's LIN- 
FOR program (see Moehler et al. 1998). 

The results are listed in Table 1 and compared to pub- 
lished values obtained from low resolution spectra. The 
formal errors of our fits arc much smaller than the sys- 
tematic errors (see below). The agreement with the results 
from low resolution spectra analysed with similar models 
(Kocn ct al. 1998) as well as from our own low resolution 
spectrum for PG 1605+072 is very encouraging. 

Four species are represented by two stages of ionization 
(He 1 and Hen, Cn and Cm, Nn and Nm, Sim and Siiv). 
Since these line ratios are very temperature sensitive at 
the temperatures in question, we alternatively can derive 
T c ff and abundances by matching these ionization equi- 
libria. Gravity is derived subsequently from the Balmer 



Table 1. Atmospheric parameters for PG 1605+072 from dif- 
ferent methods, see text 



method 


T cS [K] 


log g 


log HHa 


LTE: 








Koen et al. (1998) 


32 100+1000 


5.25+0.1 


- 


NTT, low res. 


32 300 


5.31 


-2.46 


HIRES, H+He 


31900 


5.29 


-2.54 


HIRES, He 


33 200 


5.34 


-2.42 


NLTE: 








NTT, low res. 


32 500 


5.25 


-2.46 


HIRES, H+He 


32 200 


5.21 


-2.54 


HIRES, He 


32 400 


5.21 


-2.53 


adopted 


32 300+300 


5.25+0.05 


-2.53+0.1 



lines by keeping T c s and log fixed. These two steps 
are iterated until consistency is reached. C n is represented 
by the 4267A line only, which is known to give notoriously 
too low carbon abundances. Indeed the carbon ionization 
equilibrium can not be matched at any reasonable T e g. 
The ionization equilibria of He, N and Si require T e ff to be 
higher than from the Saffer procedure, i.e. 33 200K (He), 
33900K (N) and 32 800K (Si). Table 1 lists the result for 
the He ionization equilibrium. 

This difference could be caused by NLTE effects. 
Therefore we repeated the procedure for T c ff and log ^ i 
using a grid of H-He line blanketed, metal free NLTE 
model atmospheres (Napiwotzki 1997), calculated with 
the ALI code of Werner & Drcizlcr (1999). NLTE cal- 
culations for N and Si are beyond the scope of this letter. 

Applying Saffer's procedure with the NLTE model grid 
(see Fig. 1) yields T c g almost identical to that obtained 
with the LTE grid. Evaluating the He ionization equilib- 
rium in NLTE, indeed, results in T c ff being consistent 
with that from Saffer's procedure (see Table 1). We there- 
fore conclude that the higher T c s derived above from the 
ionization equilibrium in LTE is due to NLTE effects. 

However, a systematic difference in log g persists, the 
LTE values being higher by 0.06 - 0.08 dex than the NLTE 
results (see Table 1). Since its origin is obscure, we finally 
adopted the averaged atmospheric parameters given in Ta- 
ble 1. Helium is deficient by a factor of 30 as is typical for 
sdB stars. 

4. Metal abundances 

The metal lines are sufficiently isolated to derive abun- 
dances from their equivalent widths except for the 
crowded region from 4635A to 4660A which we analyse 
by detailed spectrum synthesis. Results are listed in Ta- 
ble 2 and plotted in Figure 2. Although several O lines 
are available, it was impossible to determine the micro- 
turbulent vclocity( v t ) in the usual way, i.e. by minimizing 
the slope in a plot of the O abundances versus equivalent 
widths. Wc adopted v t ~ 5±5km/s which translates into 
small systematic abundance uncertainties of ±0.05dex for 
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most ions. The analysis is done in LTE and we therefore 
used a model ( T e ff= 33500 K, log g= 5.35) that is consis- 
tent with the He, N and Si ionization equilibria. A temper- 
ature uncertainty of A T e ff= 1000 K translates into abun- 
dance uncertainties of less than 0.1 dcx. Hence systematic 
errors are smaller for most ions than the statistical errors 
given in Table 2. 



Table 2. Metal abundances for PG 1605+072 compared to so- 
lar composition, n is the number of spectral lines per ion. 



ion 


n 


loge 


[M/H] 


Cm 


2 


7.80+0.07 


-0.83 


Nil 


5 


7.62+0.16 


-0.33 


Nm 


2 


7.66+0.24 


-0.29 


On 


25 


8.00+0.21 


-0.88 


Nen 


4 


7.57+0.19 


-0.12 


Mgn 


1 


7.16 


-0.33 


Sim 


3 


7.11+0.06 


-0.35 


Si iv 


4 


7.07+0.16 


-0.39 


Fein 


4 


7.62+0.39 


+0.30 



Carbon and oxygen are depleted by 0.8-0.9 dex with 
respect to solar composition whereas nitrogen, magnesium 
and silicon arc only slightly deficient (factor 2). Surpris- 
ingly neon and iron are solar to within error limits. This 
peculiar abundance pattern is probably due to diffusion, 
i.e. the interplay of gravitational settling and radiative 
levitation. The iron abundance is of special interest as an 
diffusive enrichment in subphotospheric layers is necessary 



to drive the pulsations. Its surface abundance is consistent 
with the diffusion calculations of Charpinet et al. (1997). 
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Fig. 2. Abundances of PG 1605+072 relative to the sun. For 
Mg an uncertainty of 0.3 dex was assumed. 



5. Projected rotational velocities 

The spectral lines of PG 1605+072 are considerably broad- 
ened, which we attribute to stellar rotation and de- 
rive vsini = 39km/s, by fitting the strongest metal 
lines. In Figure 3 we compare a section of the spec- 
trum of PG 1605+072 to that of the pulsating sdB 
star Feige48 observed with the same instrumental setup. 
However, oscillations in these objects should be asso- 
ciated with motions which could lead to line broaden- 
ing similar to the rotation effect. For a radial pulsa- 
tion with a sinusoidal velocity curve a radius change of 
AR/R=^ >3% for PG 1605+072 would be required. 
In comparison Feige48 is very sharp- lined (vsini<8km/s) 
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corresponding to AR/R<0.6% if sini=l. Although the in- 
tensity amplitudes of the pulsations in Feige 48 are consid- 
erably smaller than for PG 1605+072 we regard rotation 
as the more likely reason for the broadening of the lines of 
the latter. Time resolved spectroscopy should allow to dis- 
entangle both effects as well as to search for temperature 
changes associated with the pulsations. 

Assuming a mass of 0.5 M Q the radius of R=0.28R© 
follows from the gravity. Since sin i cannot be constrained 
the corresponding rotation period of PG 1605+072 must 
be smaller than 8.7h. PG 1605+072 displays the most com- 
plex power spectrum with more than 50 frequencies iden- 
tifiable (Kilkenny et al. 1999), 39 being bona fide normal 
pulsation frequencies. 

Usually rotation becomes manifest in the power spec- 
trum by the characteristic splitting into equidistantly 
spaced multiplct components as is observed e.g. for 
the pre-white dwarf PG 1159-035 (rotation period: 1.4d, 
Wingct et al. 1991). Such multiplct's, however, have not 
been identified for PG 1605+072. Fast rotation introduces 
higher order terms that result in unequally spaced mul- 
tiplct components. Recently, Kawaler (1999), was able to 
identify the five main peaks by considering mode trapping 
and rotational splitting. He predicted that PG 1605+072 
should be rapidly rotating (130km/s). The measured 
v sin i= 39 km/s, hence, is a nice confirmation of Kawaler's 
prediction. 

Rotation is interesting also from the point of view of 
stellar evolution. PG 1605+072 is probably already in a 
post-EHB phase of evolution (Kilkenny et al. 1999) and 
will evolve directly into a white dwarf, i.e. will shrink 
from its present radius of 0.28 R© to about 0.01 R©. Hence 
PG 1605+072 will end its life as an unusually fast rotat- 
ing white dwarf if no loss of angular momentum occurs. 
Isolated white dwarfs, however, are known to be very slow 
rotators (e.g. Heber et al. 1997). 

6. Conclusion 

Atmospheric parameters of PG 1605+072 and abundances 
of several elements have been determined at the high- 
est accuracy achieved so far for an EC14026 pulsator. 
Together with the complex pulsation pattern this makes 
PG1605+072 the most interesting object of its class. 

The high projected rotational velocity of PG 1605+072 
(v sin i=39 km/s) may be a key observation to under- 
stand its puzzling power spectrum. This measurement 
nicely confirms a recent prediction from pulsational mod- 
els (Kawaler 1999). 

Although the helium abundance is low («l/30 solar) 
as is typical for sdB stars, most metals are only mildly 
depleted (0.3-0.8 dex) whereas Ne and Fe are solar. The 
solar iron abundance is in perfect agreement with the 
predictions from the diffusion models of Charpinet et al. 
(1997). The determination of abundances of several other 
elements provides a serious test of diffusion models. 
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Fig. 3. Fit of a section of the metal line spectrum of 
PG 1605+072, (bottom, vsini=39 km/s) compared to that of 
Feige 48 (top, no rotation), another pulsating sdB star. 
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